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The on-shell imaginary part of the retarded selfenergy of massive tp 4 theory in 
1+1 dimensions is logarithmically infrared divergent. This leads to a zero in the 
spectral function, separating its usual bump into two. The twin peaks interfere in 
time-dependent correlation functions, which causes oscillating modulations on top 
of exponential-like decay, while the usual formulas for the decay rate fail. We see 
similar modulations in our numerical results for a mean field correlator, using a 
Hartree ensemble approximation. 



In our numerical simulations of 1+1 dimensional ip 4 theory using the Hartree 
ensemble approximation El we found funny modulations in a time-dependent 
correlation function. Fig. 1 shows such modulations on top of a roughly ex- 
ponential decay. The correlation function is the time average of the zero mo- 
mentum mode of the mean field, F m f(t) = ip(t)<p(0) — tp(t) (f(0), where the 
over-bar denotes a time average, X(t) — f t * dt' X(t + f')/(*2 — ti), taken after 
waiting a long time t\ for the system to be in approximate equilibrium. This 
equilibrium is approximately thermal and F m f(t) is analogous to the sym- 
metric correlation function of the quantum field theory at finite temperature, 
F(t) = (i{<p(t), <y3(0)})conn- A natural question is now, does F(t) also have 
such modulations? 

The function F(t) can be expressed in terms of the zero momentum spec- 
tral function p(p°), 

and the latter in turn in terms of the retarded selfenergy S(p°), 



p( - p '-~ (po + ie )2 + Re E (p0)]2 + [ Im E (p )]2 • U! 



The selfenergy can be calculated in perturbation theory. The one and two 
loop diagrams in the imaginary time formalism which have nontrivial energy- 
momentum dependence are shown in Fig. 2. Diagrams not shown give only 
rise to an effective temperature dependent mass, which we assume to be the 
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mass in the propagators of the diagrams in Fig. 2, after adding a counterterm 
that sets the real part of S to zero at p° — m. The one loop diagram is 
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Figure 1: Numerically computed correlation In |-F m f(t)| versus time t in units of the inverse 
temperature dependent mass m. The coupling is weak, A/m 2 = 0.11 and the tempera- 
ture T/m Kb 1.4 for the smaller volume (with significant deviations from the Bose-Einstein 
distribution) and Rb 1.6 for the larger volume (reasonable BE). 
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Figure 2: Diagrams leading to thermal damping. 



present only in the 'broken phase' (for which (<p) ^ 0; there is really only a 
symmetric phase in 1+1 dimensions, but this is due to symmetry restoration by 
nonperturbative effects which will not obliterate the one-loop damping.) The 
corresponding selfenergy has been calculated incl, for example. It only leads to 
damping for frequencies y\ > 4m 2 , which are irrelevant for the quasiparticle 
damping at Pq = m 2 . So from now on we concentrate on the two-loop diagram. 
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After analytic continuation to real time one finds that it is given by the sum of 
two terms, Si + £2 (see e.g. a). The first has an imaginary part corresponding 
to 1 «-> 3 processes requiring p 2 > 9m 2 , so it does not contribute to plasmon 
damping. The second is given by 
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where A is the coupling constant (introduced as L\ = — A<y9 4 /4), and E\ = 
v / m 2 + (p 2 +p 3 ) 2 , ^ = y/wt+tf, 1 = 2,3, n t = [exp^/T) - l]" 1 , % = 
1, 2, 3. Its imaginary part corresponds to 2 «-> 2 processes, which contribute in 
the regions near po = iro. 

Now the usual definition of the thermal plasmon damping rate (at zero 
momentum) in terms of the retarded selfenergy, 
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-ImE(m)/2m, 
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leads to a divergent answer (a collinear divergence). A natural way out of this 
difficulty may be to continue the selfenergy analytically into the lower half 

n — 77, and replace (0) by the improved 



of its second Riemann sheet, %. 
definition 

(to 
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ij) 2 + E(m - ij) = 0. (5) 

The analytic continuation of the selfenergy into the region Im p° < poses the 
puzzle how to deal with the logarithmic branch point coming from the collinear 
singularity at p° = to. However, the ambiguity is present only in the real part 



of S. For weak coupling A/m 2 <C 1 we get from (||) the equation 



9A 2 



a m/T 



167TTO 4 



( e m/T _ 1 y 



In ™ + C (T) 
7 



(G) 



The constant c has to be determined by matching a numerical evaluation of £ 
to the logarithmic singularity at p° = m. 

We evaluated £2 in @ for T = to by numerical integration with e/m = 
0.02, 0.01 and linear extrapolation e — > 0, giving c rj —0.51. For example, Eq. 
(0) now gives 7/771 = 0.061, for A/to 2 = 0.4. 

To see how well this 7 describes the decay of the correlator F(t) we eval- 
uated this function directly from ([!]) and (||). The divergence in Im£(p°) at 
p° = to leads to a zero in the spectral function p(p°). So is there a peak at 
all in p(jo 1 )? Fig. 3 shows what happens: the 'usual' peak has separated into 
two twins! Fig. 4 shows the resulting F(t). The effect of the double peak is 
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indeed an oscillating modulation on top of the roughly exponential decay. The 
decay corresponding to exp(— jt), with 7 given by (^), is also indicated in the 
plot: it does not do a good job in describing the average decay beyond the first 
interference minimum. The 'Twin Peaks' phenomenon implies that the usual 
definition of damping rate (^) is unreliable in 1+1 dimensions. 
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Figure 3: The spectral function p(p°) near p° = m = 1 corresponding to the selfenergy 
shown in Figs. 4, 5 (T = m, A = 0.4m 2 ). 
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Figure 4: Plot of In |-F(t)| versus mt for T = m, A = 0.4m 2 . The straight line represents 
exp(-Ti). 



4 



